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Abstract 

Objectives To better understand the burden of air pollution on deaths, 
we examined the effects of air pollutants on years of life lost (YLL) in 
Beijing, China. 

Design Retrospective regression analysis using daily time series. 
Setting 8 urban districts in Beijing, China. 

Participants 80 515 deaths (48 802 male, 31 713 female) recorded by 
the Beijing death classification system during 2004-08. 

Main outcome measures Associations between daily YLL and ambient 
air pollutants (particulate matter with aerodynamic diameter <2.5 u.m 
(PM 25 ), PM,„, S0 2 , and N0 2 ), after adjusting for long term trends, 
seasonality, day of the week, and weather conditions. We also examined 
mortality risk related to air pollutants. 

Results Mean concentrations of daily PM 25 , PM I0 , S0 2 and N0 2 were 
105.1 ug/m 3 , 144.6 ug/m 3 , 48.6 ug/m 3 , and 64.2 ug/m 3 , respectively. All 
air pollutants had significant effects on years of life lost when we used 
single pollutant models. An interquartile range (IQR) increase in PM 25 , 
PM 10 , S0 2 , and N0 2 was related to YLL increases of 1 5.8, 1 5.8, 1 6.2, 
and 15.1 years, respectively. The effects of air pollutants on YLL 
appeared acutely and lasted for two days (lag 0-1 ); these effects 
associated with an IQR increase in PM 25 were greater in women than 
men (11.1 (95% confidence interval 4.7 to 1 7.5) v4.7 (-2.9 to 1 2.3) YLL) 
and in people aged up to 65 years than those older than 65 years (1 2.0 
(2.9 to 21 ) 1/3.8 (-0.9 to 8.6) YLL). The mortality risk associated with 
an IQR increase in PM 25 was greater for people older than 65 years 
(2.5% (95% confidence interval 0.6% to 4.5%) increase of mortality) 
than those aged up to 65 years (0.7% (-0.8% to 2.2%)). 



Conclusions YLL provides a complementary measure for examining 
the effect of air pollutants on mortality. Increased YLL are associated 
with increased air pollution. This study highlights the need to reduce air 
pollution in Beijing, China, to protect the health of the population. 

Introduction 

The effects of air pollution on human health have recently 
attracted increasing concern in China, in part due to the 
increasing number of days with very high levels of air 
pollution. 1 2 In most Chinese cities, concentrations of PM 2 5 
(particulate matter with aerodynamic diameter <2.5 um) are 
still far above the level recommended by the World Health 
Organization' s guidelines on air quality (interim target 2 level) 
of 10 ug/m 3 (annual average) and 25 (ig/m 3 (24 h average). 3 For 
example, in 2004-08, mean daily PM 2S concentration was 105 
[ig/m 3 in Beijing. Beijing is experiencing increasing population 
density, car use, and expanded construction. It is surrounded by 
a heavy industrial region, which provides additional sources of 
air pollutants carried via air flow. Consequently, the ambient 
pollutant mixture is complex, with the potential for combined 
toxic effects from many constituents. 

Reliable estimation of the burden of air pollution on health is 
essential to support evidence based government policy in this 
important public health area. 4 5 Previous studies have examined 
the effects of air pollution on daily excess deaths or mortality 
risks using time series methods. 6 7 Those studies focused on the 
number of deaths, but did not account for age at death, apart 
from broad age stratification. We argue that using the number 
of years of life lost (YLL) provides a complementary indicator 
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to that of excess deaths, because it takes into account the life 
expectancy at death. 8 

Methods 
Data collection 

YLL data 

This study was conducted in eight districts within the urban 
area of Beijing. Mortality data on non-accidental causes were 
obtained from the death classification system at the Beijing 
Public Security Bureau, between 1 January 2004 and 31 
December 2008. These data comprised date of death, sex, and 
age. All deaths were registered residents of urban areas of 
Beijing city. 

Chinese national life tables were obtained from WHO for the 
years 2000 and 2009 (web table SI). 9 Life expectancies for 
2004-08 were averaged from the years 2000 and 2009, as data 
were unavailable for 2004-08. 

We calculated YLL for each death by matching age and sex to 
the life tables. Daily YLL were calculated by summing the YLL 
for all deaths on that day. We stratified the sums by sex and age 
group (<65 and >65 years). The web appendix shows an 
example of this calculation. 

Data on air pollution and weather conditions 

PM, 5 was monitored at the main campus of Peking University, 
located in the urban centre. 1 " " Details of the monitoring station 
are described elsewhere. 12 The monitoring station is a few 
hundred metres away from major roads and about 20 m above 
ground level. The campus is primarily residential and 
commercial without industrial sources or agricultural activities. 
Spatial variability of PM 25 mass and chemical composition is 
low across the urban area of Beijing (difference <10%). 
Additionally, average particle number and size distributions at 
this monitoring site and another regional site (50 km south of 
Peking University) were similar in the summer. 13 Therefore, the 
monitoring site provided reliable estimates of pollutant levels 
for the urban area. 12 We computed the daily average 
concentration from 24 h values. 

We obtained daily data on particulate matter less than 10 u,m 
in aerodynamic diameter (PM i0 ), sulphur dioxide (SO,), and 
nitrogen dioxide (NO,) from the Beijing municipal 
environmental monitoring centre, which had eight fixed 
monitoring sites distributed in different part of the urban area. 14 
For each monitoring site, we calculated 24 h mean 
concentrations from non-missing data if at least 18 of 24 hourly 
concentrations of PM I0 , SO-,, and N0 2 were available. 15 The 
daily mean concentrations of each air pollutant were calculated 
by averaging daily data over all monitoring sites. 

We obtained meteorological data on daily mean temperature, 
relative humidity, and air pressure from the China 
meteorological system's data sharing service. The monitoring 
station is located at Daxing district in southeast Beijing. 

Data analysis 

The daily YLL follows a normal distribution (web fig SI). We 
used daily YLL as dependent variable in a five year time series 
model, to examine its association with air pollutants. To control 
for long term trend and seasonality, we used a natural cubic 
spline with seven degrees of freedom per year for time. The day 
of week was controlled for as a categorical variable. To most 
effectively control for the potential effects of weather conditions 
on mortality, we used distributed lag non-linear models for 



temperature, relative humidity, and air pressure. A natural cubic 
spline with five degrees of freedom was used for temperature, 
relative humidity, and air pressure, and a natural cubic spline 
with four degrees of freedom for lag days (<27 days). We used 
previous experience of similar analyses in selecting the above 
parameters. 16 We validated the model fit by checking the 
residuals to ensure that seasonality and autocorrelation had been 
successfully removed. 

Studies have shown that models of single day lags might 
underestimate the effect of air pollution on mortality, 17 thus we 
used a moving average concentration over two days (lag 0-1) 
for our main analyses. 18 We also examined the associations 
using a single day lag (from lag 0 to lag 3). For each pollutant, 
we fitted models of single pollutants and multiple pollutants 
models to assess the stability of the associations. In addition, 
we stratified analyses by sex and age (<65 years and >65 years). 

To examine the linearity of the associations between air 
pollutants and YLL, we used a natural cubic spline with four 
degrees of freedom for each air pollutant (lag 0- 1 day) in single 
pollutant models. If the relations tended to be linear, we used a 
linear function; if not, we used a non-linear function with a 
natural cubic spline for air pollutants. 

To compare the standard analysis of mortality and the analysis 
of YLL, we estimated the percentage change in daily mortality 
associated with changes in air pollutants. We used the same 
independent variables as the YLL model, but with daily count 
of deaths as the dependent variable following a Poisson 
distribution. 

To check the adequacy of all models, we used an autocorrelation 
function to examine if the residuals were independent over time. 
R software was used to conduct statistical analyses." The dlnm 
package was used to perform distributed lag non-linear 
models. 20 21 

Results 

The mean concentrations of daily PM, 5 , PM I0 , SO,, and NO, 
were 105.1 u.g/m 3 , 144.6 (ig/m 3 , 48.6 (ig/m 3 , and 64.2 ug/m 3 , 
respectively (table HI). Generally, PM, 5 and PM I() had positive 
correlations with all other pollutants and weather conditions, 
while mean temperature was negatively correlated with SO, and 
N0 2 (table 2JJ). 

YLL was higher for men than women, and higher for people 
aged up to 65 years than those older than 65 years. Daily death 
counts were higher for people older than 65 years than those 
aged up to 65 years (table 1). Both YLL and death counts had 
a seasonal trend (fig 1 ij). They were higher in the cold months 
(January, February, November, and December) than the hot 
months (May, June, July, August, and September). 

The air pollutant- YLL associations tended to be linear (fig 2JJ,); 
therefore, we subsequently used a linear function for air 
pollutants. The lag pattern was similar for YLL and mortality 
risk (fig 3Jj). The effects of air pollutants appeared acutely and 
lasted for two days. Similar lag patterns were examined for 
women, men, people aged up to 65 years and those older than 
65 years (web figs S3-S6). Therefore, use of a two day moving 
average of air pollutants was sufficient to capture the short term 
effects of air pollutants on YLL and mortality risk. 

For both YLL and mortality risk, single pollutant models 
produced the highest effect estimates (table 311). For single 
pollutant models, an interquartile range increase in PM, 5 (94 
}ig/m 3 ), PM 10 (106 ng/m 3 ), SO, (49 u,g/m 3 ) and NO, (30 u,g/m 3 ) 
was related to YLL increases of 15.8, 15.8, 16.2, and 15.1 years, 
respectively. 
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The air pollutants- YLL associations differed by sex and age 
group (table 4[J). Effect estimates of PM 25 and PM I0 on YLL 
were higher in women than men, with the opposite for S0 2 and 
N0 2 . The effect estimates of air pollutants on YLL among people 
aged up to 65 years were significant and about twice those of 
people older than 65 years, although mortality risk was higher 
for older people than for those aged 65 years or younger. 

To check the robustness of our models, we performed several 
sensitivity analyses for associations between air pollution and 
YLL. All sensitivity analyses confirmed our approaches are 
valid (web appendix). 

Discussion 
Principal findings 

All air pollutants had significant effects on YLL when using 
single pollutant models, but the effect estimates decreased when 
multiple pollutants were included in the models. Effects 
appeared acutely and lasted for two days (lag 0-1). The effects 
of air pollutants on YLL were higher in people aged up to 65 
years than those older than 65 years, whereas these results were 
opposite for death counts. 

Interpreting the findings 

One crucial finding was that the estimated effects of air pollution 
were greater on the younger group of people than older group 
for YLL. The potential reason is that the measurement of YLL 
takes into account those conditions afflicting young people or 
children. Giving the same weight to deaths occurring at different 
ages could distort policy priorities and resource allocation. 22 
Most studies report that mortality risks related to air pollution 
are greater for older people than younger people. 23 24 Our study 
suggests that focusing on death counts only could underestimate 
the burden of air pollution on young people. 

Many studies have reported that the effects of air pollution on 
women are higher than men. Both biological and non-biological 
factors are associated with this difference. Women have smaller 
lungs and airway diameters. These might increase airway 
reactivity and exacerbate particulate deposition. 25 Women and 
men have different socioeconomic status and stress 
experiences. 26 27 Women also tend to spend more time outside 
because they are less likely to get full time jobs. 28 

Our results show that the harvesting effects existed only in older 
people and women for both YLL and mortality risk. Previous 
studies have confirmed this pattern. 29 30 The public health 
significance of air pollution becomes smaller, if there is 
mortality displacement. However, we did not find harvesting 
effects in men or younger people. 

All pollutants have considerable effects on YLL and mortality 
risk in Bejing, which implies that stricter standards should be 
in place for all air pollutants. When we used the two or three 
pollutant models, effects of air pollutants were reduced. These 
findings were consistent with previous studies, 14 18 24 and may 
be caused by co-linearity between air pollutants, in turn caused 
by commonality of sources or photochemical interactions. 

Meaning of the study 

The findings strongly support the need for authorities to reduce 
air pollution in Beijing. Although the general population is 
increasingly concerning with air pollution, information is still 
limited. Real time data on PM 25 have been available since 2012, 
after a public outcry about thick smog in China. Emergency 
response measures have been implemented on high smog days, 
including reducing industrial emissions, removing some 



government vehicles from the road, and halting outdoor 
activities for school children. The effective reduction in 
emissions during the 2008 Beijing Olympics demonstrates that 
it is possible to reduce air pollution in China. 

Strengths and limitations of the study 

Our study had some limitations. We used ambient pollutant 
concentrations as surrogates of individual exposure, which could 
result in measurement error. Because data were from only one 
city, it is difficult to generalise results to other cities. We did 
not control for smoking or for prevalence of chronic obstructive 
pulmonary disease related to smoking, because this information 
was unavailable. However, we assume that our results would 
have been little affected because the effects are short term and 
the time series method controls for long term and fixed term 
factors, such as smoking and obesity. The distribution of such 
fixed factors does not vary from day to day, and thus they are 
not associated with air pollution levels. 31 

This study also had some strengths. It examined the burden of 
air pollutants on YLL in China. Compared with mortality risk 
that weighs all deaths equally, YLL is a more informative 
indicator for quantifying premature deaths. Our findings are 
important in developing public policy, determining appropriate 
interventions to manage risk, and promoting capacity building 
for local responses to air pollution. 

Conclusions 

This study highlights the effects of exposure to air pollution on 
YLL and mortality risks in Beijing, China. Our findings support 
the need to reduce the high levels of air pollution in Beijing, 
China. YLL can be used as a complementary indicator for 
assessing the effect of air pollutants on mortality. 
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What is already known on this topic 

Air pollution increases the risk of mortality, and is a serious problem in Beijing, China 

However, no study has examined the burden of air pollution in terms of years of life lost, which combines the counts of deaths with life 
expectancy 

What this study adds 

Between 2004 and 2008, all air pollutants (PM 25 , PM 10 , S0 3 , and NOJ had significant effects on years of life lost in Beijing; the effects 
of air pollutants on years of life lost appeared acutely and only lasted for two days (lag 0-1) 

People aged up to 65 years were more affected by air pollutants than those older than 65 years for years of life lost, while the mortality 
risk was higher for those older than 65 years than those aged up to 65 years 

Years of life lost is more informative for quantifying premature deaths than mortality risk, which weighs all deaths equally 



aspects of the study have been omitted; and that any discrepancies 
from the study as planned (and, if relevant, registered) have been 
explained. 
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Tables 



| Levels of dally air pollutants, mean temperature, relative humidity, mean air pressure, and YLL In Beijing, China, 2004-08 





Minimum 25% quartile Median 75% quartile Maximum 


Mean 


Standard deviation Interquartile range 


PM 25 (ug/m 3 ) 


0.7 


45.3 


85.3 


139.2 


517.7 


105.1 


80.9 


94 


PM t0 (ug/m 3 ) 


10.0 


78.3 


128.0 


184.0 


600.0 


144.6 


91.3 


106 


S0 2 (ug/m 3 ) 


5.0 


15.0 


29.0 


63.8 


293.0 


48.6 


49.1 


49 


N0 2 (ug/m 3 ) 


14.0 


47.1 


60.9 


76.8 


214.4 


64.2 


25.7 


30 


Mean temperature (°C) 


-10.1 


3.5 


14.9 


23.6 


32.1 


13.6 


10.9 


20 


Relative humidity (%) 


8.0 


34.0 


52.0 


68.0 


97.0 


51.5 


20.3 


34 


Air pressure (Pa) 


987.8 


1004.0 


1012.0 


1021.0 


1043.0 


1013.0 


10.2 


17 


YLL (years) 


Total 


299.4 


593.8 


695.9 


819.1 


1504.0 


709.2 


160.5 


225 


Women 


52.6 


224.1 


282.5 


349.5 


601.9 


289.2 


91.8 


125 


Men 


117.1 


342.7 


412.4 


491.9 


975.4 


420.0 


110.7 


149 


Age <65 years 


105.9 


323.5 


407.3 


498.2 


954.7 


416.4 


127.2 


175 


Age >65 years 


75.6 


241.2 


286.6 


340.1 


605.5 


292.8 


73.4 


99 


Daily death counts (No of deaths) 


Total 


19 


38 


43 


50 


86 


44.1 


9.1 


12 


Women 


5 


14 


17 


21 


36 


17.4 


4.9 


7 


Men 


9 


22 


27 


31 


53 


26.7 


6.4 


9 


Age <65 years 


4 


12 


15 


18 


34 


15.2 


4.5 


6 


Age >65 years 


7 


24 


28 


33 


56 


28.8 


7.1 


9 
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Table 2| Spearman correlation between air pollutants and weather conditions in Beijing, China, during 2004-08 





PM 10 


so 2 


N0 2 


Mean temperature Relative humidity Air pressure 


PM, 5 


0.67* 


0.32* 


0.61* 


0.07 


0.39* 


-0.13* 


PM,„ 




0.44t 


0.67* 


0.05 


0.21* 


-0.13* 


S0 2 






0.60* 


-0.71* 


-0.26* 


0.58* 


N0 2 








-0.16* 


0.18* 


0.12* 


Mean temperature 










0.36* 


0.87* 


Relative humidity 












-0.35* 


*P<0.05. 
tP<0.01. 
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Table | Association between interquartile range increase in air pollutants (lag 0-1 day) and YLL and increase in deaths for non-accidental 
deaths using single, two, and three pollutant models during 2004-08 



Pollutant^ and model 


YLL (years) 


Increase in deaths (%) 




Single model 


15.8 (5.3 to 26.3)t 


1.3(0.1 to 2.6)* 


+S0 2 


9.7 (-3.0 to 22.4) 


0.5 (-1.0 to 2.1) 


+N0 2 


7.8 (-6.3 to 21.9) 


0.3 (-1.4 to 2.0) 


+S0 2 +N0 2 


6.9 (-7.4 to 21.1) 


0.2 (-1.5 to 1.9) 




Single model 


15.8 (6.1 to25.5)t 


1.7 (0.6 to 2.9)f 


+S0 2 


12.6 (0.7 to 24.4)* 


1.3 (-0.1 to 2.7) 


+N0 2 


11.9 (-2.1 to 25.9) 


1 .3 (-0.4 to 3.0) 


+S0 2 +N0 2 


11.6 (-2.5 to 25.6) 


1.3 (-0.4 to 3.0) 


S02 


Single model 


16.2 (4.1 to28.4)t 


1 .8 (0.4 to 3.2)t 


+ PM,5 


14.2 (-2.4 to 30.8) 


1 .7 (-0.2 to 3.6) 


+PM t0 


7.1 (-7.9 to 22.0) 


0.8 (-0.9 to 2.6) 


+N0 2 


6.8 (-11. 7 to 25.3) 


1.0 (-1.2 to 3.1) 


+PM 25 +N0 2 


8.9 (-1 1 .4 to 29.3) 


1.1 (-1.2 to 3.5) 


+PM 10 +NO 2 


5.6 (-13.0 to 24.2) 


0.9 (-1.3 to 3.0) 


N02 


Single model 


15.1 (4.7 to 25.6)t 


1.6 (0.4 to 2.8)t 




13.2 (-2.3 to 28.7) 


1 .6 (-0.2 to 3.5) 


+PM 10 


5.8 (-9.3 to 21.0) 


0.5 (-1.3 to 2.3) 


+S0 2 


10.8 (-5.1 to 26.6) 


0.9 (-1.0 to 2.8) 


+PM 25 +S0 2 


8.4 (-10.6 to 27.3) 


1.0 (-1.3 to 3.3) 


+PM, 0 +SO 2 


2.5 (-16.3 to 21 .3) 


0.0 (-2.3 to 2.2) 



Data are mean (95% confidence interval) and are controlled for seasonality, day of the week, temperature, relative humidity, and air pressure. 

*P<0.05. 

tP<0.01. 

^ Interquartile ranges were 94 ug/m 3 for PM 25 , 106 ug/m 3 for PM I0 , 49 ug/m 3 for S0 2 , and 30 ug/m 3 for NO z . 
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Table | Association between interquartile range increase in air pollutants (lag 0-1 day) and YLL and increase in deaths for non-accidental 
deaths using single pollutant models during 2004-08, according to sex and age 





Sex 




Age 




Pollutant* 


Female 


Male 


<65 years 


>65 years 


YLL (years) 




11.1 (4.7 to 17.5)t 


4.7 (-2.9 to 12.3) 


12.0 (2.9 to 21 )t 


3.8 (-0.9 to 8.6) 


PM 10 


9.3 (3.3 to 15.2)f 


6.5 (-0.5 to 13.5) 


10.3 (2to18.6)t 


5.5(1.1 to 9.9)* 


S0 2 


5.6 (-1.9 to 13.1) 


10.6 (1.8 to 19.4)* 


10.8 (0.3 to 21.3) 


5.4 (-0.1 to 10.9) 


NO, 


6.7 (0.3 to 13.1)* 


8.4 (0.8 to 16.0)* 


10.1 (1.1 to 19.1)* 


5.0 (0.3 to 9.8) 


Increase of deaths (%) 




2.2 (0.4 to 4.1)* 


0.8 (-0.8 to 2.3) 


0.7 (-0.8 to 2.2) 


2.5 (0.6 to 4.5)* 


PM 10 


2.5 (0.8 to 4.2)t 


1.2 (-0.2 to 2.6) 


1.3 (-0.1 to 2.7) 


2.5 (0.6 to 4.4)* 


S0 2 


1.9 (-0.2 to 4) 


1.7 (-0.1 to 3.4) 


1 .3 (-0.4 to 2.9) 


2.8 (0.5 to 5.1)* 


NO, 


1.9 (0.1 to 3.7)* 


1.4 (-0.2 to 2.9) 


1 .2 (-0.3 to 2.6) 


2.4 (0.4 to 4.4)* 



Data are mean (95% confidence interval) and are controlled for seasonality, day of the week, temperature, relative humidity, and air pressure. 

*P<0.05. 

tP<0.01. 

^ Interquartile ranges were 94 ug/m 3 for PM 25 , 106 ug/m 3 for PM, 0 , 49 ug/m 3 for S0 2 , and 30 ug/m 3 for NO z . 
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Figures 

Total Female Male <65 years >65 years 




Fig 1 Boxplots of monthly YLL and death counts in Beijing, China, during 2004-08, according to sex and age 
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Fig 2 Association between air pollutants (lag 0-1 day) and YLL in Beijing China, during 2004-08. A natural cubic spline with 
four degrees of freedom for air pollutants was included in the single pollutant models, while controlling for seasonality, day 
of the week, temperature, relative humidity, and air pressure 
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Fig 3 Association between increased interquartile range in air pollutants and YLL (top) and percentage increase of deaths 
(bottom) for non-accidental deaths using single pollutant models at different lag days, during 2004-08. Results were controlled 
for seasonality, day of the week, temperature, relative humidity, and air pressure. Interquartile ranges were 94 ug/m 3 for 
PM 25 , 1 06 ug/m 3 for PM 10 , 49 ug/m 3 for S0 2 , and 30 ug/m 3 for N0 2 
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